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Noticeable deviations from the prediction of the fiducial LCDM cosmology are found in the angular 
power spectrum of the CMB. Besides large-angle anomalies, the WMAP 1*' year data revealed a 
dip in the power spectrum at I ~ 200, which seemed to disappear in the 3'^'* year and subsequent 
angular power spectra. Using the WMAP single 1°*, 3'^'', and S**^ year data as well as the total 5 
year coadded data, we study the intensity and spatial distribution of this feature in order to unveil 
its origin and its implications for the cosmological parameters. We show that in all WMAP data 
releases there is a substantial suppression of the first Doppler peak in a region near the north ecliptic 
pole. 



I. INTRODUCTION 

The dynamics of both individual galaxies and clusters 
of galaxies, the Hubble diagram for distant Type la su- 
pernovae and the pattern of temperature anisotropics in 
the Cosmic Microwave Background (CMB) broadly sup- 
port the standard flat A-dominated cosmological model 
with a nearly scale- invariant spectrum of adiabatic Gaus- 
sian primordial fluctuations, such as might be generated 
at the end of an inflationary epoch, as an accurate de- 
scription of our universe. 

In particular, this simple model provides an ac- 
ceptable flt to the Wilkinson Microwave Anisotropy 
Probe (WMAP) observations of the CMB temperature 
anisotropics. The pattern of acoustic peaks and troughs 
predicted by the model to be imprinted on the angu- 
lar power spectrum by the primordial fluctuations in the 
inflationary field and its subsequent evolution has been 
detected and characterized to good precision [TH3]. A 
global fit to the data has enabled cosmologists to deter- 
mine the handful of parameters of the inflationary ACDM 
model, which does a remarkably good job at describing 
the spectrum. 

The WMAP data, however, also exhibits unexpected 
anomalies that have sparked considerable attention |4| 
[S]. These include excesses or deficiencies of power in 
at least three bins - £ ~ 22, 44 and 200 - violations 
of statistical isotropy at £ < 6 [Ml], and a severe lack 
of large angle correlations requiring the low-£ C^ to not 
be independent [5J I19j . Hemispheric asymmetries in the 
temperature power spectrum extending over a wide range 
of angular scales as well as hemispheric dependent non- 
gaussian signatures were reported in [9Ul2j. and a non- 
Gaussian cold spot in the southern galactic hemisphere 
was found using wavelet transforms |14l I15j . Evidence 
for primordial non-Gaussianity of the local type in the 
temperature spectrum in the third year data has been 
claimed in [TC| (However, it does not seem to be present 



when later data are included in the analysis p71ll8j ). 

If these anomalies are of primordial origin, rather than 
caused by unaccounted foreground or secondary effects, 
there could be profound implications for our understand- 
ing of the early Universe. For example, while the simplest 
single field infiationary scenarios predict negligible devi- 
ations from Gaussianity, larger deviations can be accom- 
modated in extended particle physics models j21j . The 
lack of power in some low-£ multipoles, though not the 
multipole alignments and not the lack of large-angle cor- 
relations, could be explained by the presence of anticorre- 
lated isocurvature perturbations such as those expected 
in the curvaton model [^ [25| . Topological defects could 
be at the origin of the cold spot [21] . 

Most of these anomalies were first noted in the first 
year WMAP data release and persisted into the subse- 
quent data releases [ini HO]- The dip in the first peak 
found in the first year WMAP data release is believed to 
be a notable exception - it is no longer present when the 
newly adopted noise weighting scheme is used to extract 
the best flt power spectrum. This anomaly in the flrst 
year angular power spectrum was, thus, attributed to a 
noise fluctuation [3]. 

Surprisingly, the dip in the power spectrum around 
£ ~ 200 disappears when data from the ecliptic poles 
is not included in the analysis of the 1st year WMAP 
data regardless of the weighting scheme (see Fig. 7 
in P). Archeops, the sole previous ground-based ex- 
periment that observed the region of the sky around the 
north ecliptic pole, showed a dip around the first peak. 
All other similar experiments showed no such dip and 
had access only to southern ecliptic regions [3S]. This 
suggests that the differences might not simply be due to a 
noise fluctuation, but might rather reflect a real anomaly. 

In this paper we take a closer look at the degree-scale 
CMB temperature anisotropics. We estimate the pri- 
mordial power spectrum using data from different re- 
gions of the sky according to the two weighting schemes 
that have been used by the WMAP collaboration. Our 



findings show that there is a substantial suppression of 
the first Doppler peak in a region near the north ecliptic 
pole. This reduction in the temperature power spectrum 
around the first acoustic peak occurs regardless of the 
weighting scheme used, although the detailed shape (i.e. 
the dip like feature) is only present when the analysis 
follows the strategy used by WMAP in their first year 
release and might also depend on the particular binning 
used to present the data. 

The rest of the paper is organized as follows. In sec- 
tion [nj we review the method used for estimating the 
power spectrum, Cg, from observations. Section |lIT] de- 
scribes the datasets that we used, together with the pro- 
cessing steps taken to minimize the influence of non- 
cosmological foregrounds. In section |IV| we discuss the 
method used to determine the statistical significance of 
our analysis of the data. Next, in section |VJ before our 
conclusions, we present the results of our analysis and 
briefly comment on implications for the physics underly- 
ing the CMB. 



II. ESTIMATION OF POWER SPECTRA 

Two main approaches have been used to estimate the 
power spectrum from a CMB data map, including sources 
of errors and partial sky coverage, namely maximum- 
likelihood estimators and pseudo-C^ (PCL) methods. 
The WMAP team used the latter strategy for the analysis 
of their first year data, but changed to a hybrid estimator 
when producing the third year results, in agreement with 
the discussion in |26| . The differences, however, concern 
only the lowest multipoles {i < 30), which are not the 
object of our study. Hence, we base our results on the 
pseudo-C^ technique |27) . which we now briefly review. 

The usual decomposition of a sky map AT(n) in spher- 
ical harmonics can be generalised to the case of partial 
sky coverage by introducing a position-dependent weight, 
W{n), which is set to zero in the regions of the sky that 
are to be excluded from the analysis: 
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Here the map has been discretized in pixels sub-tending 
a solid angle f2p, and the index i refers to the particu- 
lar differential assembly (DA) under scrutiny (see [I] for 
details on the datasets). 

Even though the pseudo-power spectrum, Ci, derived 
from weighted maps M, 
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differs from the full-sky angular spectrum, C| ^, their 
ensemble averages are related by a mode coupling matrix, 
Gw , which can be inverted to obtain an estimator of the 
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The coupling matrix depends only on the geometry of 
the weight function, VF(n). An expression for VK(n) can 
be found in j27j . 

In a multichannel experiment like WMAP, the noise 
between two different DAs is assumed to be uncorrelated. 
With a straightforward generalization of Eq. (pi) for the 
analysis of two DA maps at a time. 
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we obtain an estimate of the true power spectrum that 
is not biased by noise. 

The weight function, VF(n), is zero where the fore- 
ground emission cannot be reliably eliminated. The main 
sources of non-cosmological origin (free-free, dust and 
synchrotron emission) are subtracted from each DA map 
by fitting templates and taking advantage of the differ- 
ent frequency dependence of each signal. The intensity 
of the foregrounds in the plane of the Galaxy, however, 
disallows a proper separation of the cosmological compo- 
nent in this region, which should be excluded from our 
analysis. We will set W to zero in the region specified by 
the KQ85 five year mask, an area of roughly 15% of the 
sky. 

It is also advantageous to choose a different form of 
VF(n) in those regions where it is non-zero, depend- 
ing on the multipole angular scale involved. Maximum- 
likelihood estimators, leading to smaller error bars, 
weight the data by the inverse covariance matrix, C^^ — 
{S + N)~^, where S and N are the contributions from 
the signal and the instrument noise respectively. We can 
mimic this optimal sampling for the PCL estimator by 
using a unit uniform weight in the signal dominated, 
C^^ ~ S^^, low £ regime. An inverse noise weight, 
WP sa N^i,^, is more suitable for the noise dominated 
high multipoles [26]. In the first-year WMAP analysis, 
the two regimes were taken to be £ < 200 and £ > 450, 
while a transitional weighting interpolation was used for 
the intermediate multipoles |^ . For the third year anal- 
ysis, the WMAP team used a sharp transition between 
the uniform weight and inverse noise weight at £ = 500 
(the transition was made at £ = 600 for the seventh year 
analysis, although this does not affect our discussion). 
The fact that the dip-like feature in the first acoustic 
peak is not evident in the WMAP third year plot, lead 
the WMAP collaboration to conclude that it was largely 
due to the sharp weight transition at £ < 200. We will see 
in the following that, in addition to this effect, the CMB 
power spectrum at these scales seems to vary markedly 
among different regions of the sky. 



III. DATA SETS 

The WMAP team obtained its best estimate power 
spectrum by optimally weighting all the possible pairings 
of Q, V and W DAs, although the Q band was dropped in 
the third year estimate, since it is the most prone to fore- 
ground and diffuse source contamination. Our aim here 
is not to obtain the best power spectrum estimate, but to 
study how it changes when looking at different parts of 
the sky or when using different weighting schemes. We 
thus obtain our power spectrum estimate by applying 
Eq. (|3| to the foreground reduced V and W frequency 
band maps for the single first, third, and fifth year maps 
as well as the complete fifth year coadded map. We per- 
form our analysis using the data processing method cor- 
responding to each data release year as well as the ap- 
plying the first and third year weighting schemes [TJ |3] 
to all datasets. 

Oddly enough, some of the anomalies found at large 
angular scales correlate with the ecliptic [7|, and the 
WMAP team first year data gives a different angular 
power spectrum around (. ~ 200 in the ecliptic plane 
than at the poles [T]. A random search in all directions 
could potentially reveal a more significant region, how- 
ever when such an analysis was performed for the 1st 
year data the results also pointed to the ecliptic poles, 
as Eriksen, et al. found significantly low power in the 
northern hemisphere near the ecliptic pole for the range 
£ = 2 — 40 1^. We, thus, choose to examine two a pri- 
ori selected 30 degree diameter spherical caps centered 
on the ecliptic poles. A later analysis (also on the first 
year data) that computed the directional dependence of 
the temperature power spectrum for an € range that in- 
cluded the first peak was done in |13j . A comparison of 
those results with the results of our analysis is left for 
Section El 

Estimating the real power spectrum from observations 
of a limited region in the sky, introduces errors that in- 
crease as we shrink the patch. To minimize this effect 
we use complementary masks. For instance, to check the 
influence of the northern ecliptic pole region, we just re- 
move the cap around the pole and keep the rest of the 
sky. 



IV. ESTIMATION OF STATISTICAL 
SIGNIFICANCE 

Masking a full sky map introduces correlations between 
the different multipoles that limit the precision of the es- 
timated underlying power spectrum. Because our anal- 
ysis of the data includes masking large regions near the 
ecliptic poles, we must determine whether our results are 
statistically significant or just a result of masking the 
map. 

To do this, we generate IC* random realizations of the 
sky map from the best fit ACDM power spectrum. For 
each data release year and frequency band we wish to 



simulate we apply the corresponding beam transfer func- 
tion, supplied by the WMAP team, to the extracted aem- 
Next, we generate two identical maps of the sky and 
add uncorrelated noise to the pixels of each, given by the 
function 



Noise — aai{p), 



where 
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is the noise for each pixel p in the ith frequency band map 
with i = W or V. Here a is a randomly generated number 
with a Gaussian distribution, cto.j the mean rms noise 
used to calculate the pixel noise for the zth band provided 
by [T] , and Nobs (p) is the number of observations recorded 
for pixel p that is taken from the corresponding data 
release year and frequency band WMAP data map we 
wish to simulate. This mimics the noise present in the V- 
and W-band data maps that, being independent, drops 
from the cross-correlation in Eq. Q. 

Each map is then masked; first only the WMAP KQ85 
mask is used, then the analysis is repeated with an addi- 
tional region masked out (these regions include the north, 
south ecliptic poles and/or the ecliptic plane). We sub- 
tract the monopole and dipole, which are induced during 
masking, and calculate the cross-power spectrum using 
Eq. (HI). Finally, we remove the window function and 
calculate the true CiS from the pseudo-C^s using Eq. (|3|. 

We would like to quantify the relative difference be- 
tween a map masked with only the KQ85 mask and the 
same map masked with KQ85 -f- some additional region. 
We have thus defined a peak statistic 

Eima^c f ^iKQSb+add'l mask ^KQS5\ 
_ £=£„,„ \^£ ~ '-'e J 

peak — «^^^ ( f^KQS5+add'l mask ^KQS5\ ' ^ ' 

Here imin and imax are the extrema of a range of £ that 
covers the first peak of the power spectrum. The exact 
^mjTi=200 and i!maa;=250 values were chosen such that the 

value of imax {imax + ^)Ce^^^ and imin{imin + ^)Cl^.^ 

are 90% of the WMAP five year power spectrum peak 
value at its maximum. 



V. RESULTS AND DISCUSSION 

The observed power spectrum in the first year WMAP 
data for scales i ^ 200 was ^10% smaller when using 
only the data from the ecliptic poles than what the best 
fit ACDM model would indicate. If this dip originated 
due to the lack of power entirely in a particular region in 
the sky, we can then estimate the diameter of the circular 
patch involved to be ~ 30°. 

As shown in Fig. Ill the power around £ ^ 200 in- 
creases if data in the region of the north ecliptic pole is 
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{a)First year data. 



(a)First year data. 




6500 










/' - ^' '■ . 




- 


„ 6000 








,y 


^"""""^^^ij-T ^ 


\ 




_3, 








5- 5500 
+ 




rf/ 










- 


5000 












, \\ 





160 180 200 220 240 260 280 300 



160 180 200 220 240 260 280 300 



(b)Third year data. 



(b)Third year data. 
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(c)Fifth year data. 



(c)Fifth year data. 



FIG. 1: Power spectra for each data release year using the 
first year weighting schemes. We show the power spectra 
extracted from the fuU sky (blue solid), full sky without the 
north (red dashed) and south (black dot dashed) ecliptic caps, 
the ecliptic plane only (green dotted), and the ecliptic poles 
only (yellow circles). The galaxy is always masked out. 



FIG. 2: Power spectra for each data release year using the 
third year weighting schemes. We show the power spectra 
extracted from the full sky (blue solid), full sky without the 
north (red dashed) and south (black dot dashed) ecliptic caps, 
the ecliptic plane only (green dotted), and the ecliptic poles 
only (yellow circles). The galaxy is always masked out. 
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(a)First year weighting. 




160 180 200 220 240 260 280 



(b)Tliird year weighting. 



FIG. 3: Power spectra for the full 5-year coadded data using 
both the first and third year weighting scheme. We show 
the power spectra extracted from the full sky (blue solid), 
full sky without the north (red dashed) and south (black dot 
dashed) ecliptic caps, the ecliptic plane only (green dotted), 
and the ecliptic poles only (yellow circles). The galaxy is 
always masked out. 



masked out when estimating the spectrum. The blue solid 
line was obtained with the KQ85 mask only, and for the 
black dot dashed line we left out, in addition, the south- 
ern ecliptic pole. When leaving out the northern ecliptic 
pole (red dashed) or keeping just the ecliptic plane (green 
dotted) the value of Cf^200 is roughly 10% higher. 

This effect is seen for all WMAP data analyzed here 
regardless of the chosen weighting scheme as shown in 
Fig. [2] for the third year weighting scheme analysis and, 
more importantly, in Fig. [3] where we have used the full 
5 year coadded data. Although, the dip-like shape in the 



full sky minus the galaxy only appears when the first year 
data and weighting scheme is used. The 5 year coadded 
data is the most compelling part of the analysis, since 
there are far more observations per pixel and therefore 
less noise is in the resulting maps. 

The significance of this effect can only be assessed with 
a detailed likelihood analysis. To this end, we have cal- 
culated the p-value for the peak statistic, given by Eq. 7 
of each data release and weighting scheme. Figs. |4j 5 
& |6] show the probability densities of peak statistic val- 
ues found for the simulations masked with the KQ85 and 
north ecliptic-|-KQ85 masks for the first and third year 
weighting schemes respectively. The peak statistics and 
corresponding p-values (shown as percentages) for each 
data set are listed in Table |l] These percentages are cal- 
culated by integrating under the curve in Figs. |4] [5j &[6] 
to the left of the specific value of the peak statistic found 
for the data. The significance is above ~95% regardless 
of the data set and weighting scheme. 

The process of foreground cleaning was overhauled in 
the third year WMAP data release, and could account for 
the few percent drop in the significance value from the 
first to third year analysis. Still, the reduction in power 
due to northern ecliptic pole data is still significant (by 
about 2a) in the later data releases. 

As for possible systematic effects, it should be noted 
that the WMAP surveys the ecliptic plane more sparsely 
than the poles, and it has been suggested that beam ellip- 
ticity could bias the result. However, Wehus et al. [29] 
performed an extensive analysis of the effects of asym- 
metric beam patterns on the power spectrum, and found 
the change to be much smaller than what is shown in 
Figs.[lJ[2]&[3]. 

Hansen, et. al |13j performed an analysis on the first 
year WMAP data by varying the location of 9.5° and 
19° disks around the sky and computing the power spec- 
trum. They found that 3 disks resulted in a 2a result, 
however this happened in 10% of their simulations which 
lead them to conclude that the effect was not signifi- 
cant. One important difference between the Hansen, et. 
al analysis and our own is that they used only the map 
data from inside the disk and consequently had use large 
£ bins to reconstruct the power spectrum. Our method, 
which uses complementary masks - analysis is on the 
data outside of the disc rather than inside, allows us to 
use a much smaller range of £ for each bin. Therefore, 
we are able to resolve the first peak in the power spec- 
trum with a higher degree of accuracy (we have 6 points 
for 150 < £ < 250 whereas they have 2). We believe 
this could account for the difference in the determined 
significance of this effect. 



VI. CONCLUSION 

In the first year WMAP TT power spectrum, the 
binned C'l value at the first peak was significantly lower 
than the best-fit ACDM value [T]. This was widely re- 
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FIG. 4: Probability density of the peak statistics calculataed 
using first 4(a) third 4(b) and fifth 4(c) year simulated data 



FIG. 5: Probability density of the peak statistics calculataed 
using first 5(a) third 5(b) and fifth 5(c) year simulated data 



for the first year weighting scheme. The y-axis shows the 
normalized number of Monte Carlo skies that have each par- 
ticular value of S. 



for the third year weighting scheme. The y-axis shows the 
normalized number of Monte Carlo skies that have each par- 
ticular value of S. 
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FIG. 6: Probability density of the peak statistics calculataed 



using the five year coadded simulated data with the first 6(a) 



third 6(b) year weighting schemes. The y-axis shows the nor 
malized number of Monte Carlo skies that have each particu 
lar value of S. 



TABLE I: The first columns under the weight scheme head- 
ings list the peak statistic calculated for each data release. 
The second columns list the corresponding significance. 
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ferred to as the "dip in the first peak." In the third year 
data release the so-called dip disappeared. This was at- 
tributed to the use of an improved weighting scheme [3] . 
However, there was also shown to be a decrease in power 
at the first peak from the part of the sky near the eclip- 
tic poles compared to the ecliptic plane in the first year 
data [T] . Indeed, using the improved (third year) weight- 
ing scheme, we have found that the ecliptic polar power 
at the first peak is significantly lower than the full-sky 
power in the first, third, fifth and coadded 5 year data 
releases. (We have not yet analyzed the seventh year data 
release, although we do not expect much difference.) 

By separately masking out the north ecliptic polar re- 
gion, the south ecliptic polar region and the ecliptic plane 
we have demonstrated that this difference can be ascribed 
to a large reduction in peak power in the north ecliptic 
polar region. We have shown that there is a significant 
increase in power when we mask out the north ecliptic 
pole for all weighting schemes and data releases. There is 
no such change in power when the south ecliptic pole is 
masked. Additionally, when using the third year weight- 
ing scheme there is still a sizeable reduction in power 
when only the data from the ecliptic poles is used com- 
pared to when the full sky is analyzed. 

We have characterized the absence of power at the first 
doppler peak in the temperature power spectrum in the 
north ecliptic polar region using the peak statistic defined 
in Eq. [7] Using either the first or third year weighting 
scheme and across all data releases the significance is at 
least 95% based on 10'* best-fit ACDM realizations and 
the analysis of one V x W cross-band power spectrum. 
The most trusted result should arguably be the signif- 
icance of the coadded map with the third year weight 
scheme, as it provides analysis with the least noise and 
data processing method closest to what is currently used. 
Even when analyzing the most optimal data set for char- 
acterizing this effect, our result is still « 2a significant. 

Indeed, the first peak may seem like a specific choice of 
range in i to study - we would like to extend this analysis 
in the future to other "special" regions in £-space, namely 
to the second and third doppler peaks. However current 
experiments are limited in their measurements around 
these peaks, so an accurate characterization of this effect 
is not possible. Future CMB experiments should provide 
better data for further analysis. 

Archeops's observation of a dip in the first peak from 
measurements including the north ecliptic polar region of 
the sky, albeit at low statistical significance, suggests that 
this is not a WMAP systematic. Planck should therefore 
confirm the existence of a dip in the first peak power near 
the north ecliptic pole. 

If these variations are of cosmological origin, they can- 
not be explained within the standard ACDM model, since 
it predicts an isotropic and Gaussian spectrum. How- 
ever, the association with the north ecliptic pole argues 
against a cosmological origin and in favor of Solar Sys- 
tem physics. Possible signatures include a deviation in 
the CMB spectrum in this region of the sky and an as- 



sociated anomaly in polarization maps. 
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